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Abstract 

Fluids often display dissipative properties. We explore dissipation in the form of bulk 
viscosity in the cold dark matter fluid. We constrain this model using current data from 
supernovae, baryon acoustic oscillations and the cosmic microwave background. Considering 
the isotropic and homogeneous background only, viscous dark matter is allowed to have a bulk 
, viscosity < 10 7 Pa-s, also consistent with the expected integrated Sachs- Wolfe effect (which 

plagues some models with bulk viscosity). We further investigate the small-scale formation of 
viscous dark matter halos, which turns out to place significantly stronger constraints on the 
dark matter viscosity. The existence of dwarf galaxies is guaranteed only for much smaller 
values of the dark matter viscosity, < 10~ 3 Pa-s. 

PACS numbers: 98.80-k, 95.36+x 

q . 1 Introduction 

Modern cosmology assumes the existence of dark energy in the form of a cosmological constant A 
^ ', to dominate the cosmic energy budget today. Dark energy is believed to cause the current accel- 

ii erated expansion of the Universe. The matter components of the Universe are baryons, photons, 

neutrinos and cold dark matter (CDM), which are modelled as independent and dissipationless 
. fluids. 

This concordance model assumes from the very beginning the cosmological principle and 
relies on the validity of the first law of thermodynamics. This means that directional dissipative 
processes, i.e. heat conduction, diffusion and shear stresses, are excluded due to the high degree 
of isotropy imposed to the Universe. However, bulk viscosity, which is induced by a divergence 
of the velocity field, is the unique viscous mechanism allowed for cosmic fluids even in a ho- 
mogeneous and isotropic background. Hence, the effective pressure of any cosmic fluid can be 
written as the sum of the kinetic pressure, p k , and the bulk viscous pressure, p v , such as (c = 1) 



o 
o 



Pes =Pk + Pv = w k e - 3H£, (1) 



where the coefficient of bulk viscosity £ is a non-negative quantity £ > 0, due to the second law 
of thermodynamics It is obvious that the bulk viscosity reduces the effect of the kinetic 
pressure and can even make the effective pressure negative (see e.g. [2]). Assuming ([1]), the 
continuity equation of a bulk viscous fluid reads 

e + 3H(e + p k )-9H 2 £ = 0. (2) 

The purpose of this paper is to investigate whether observations allow the presence of bulk 
viscosity in the dark matter fluid. 

Clusters of galaxies revealed the existence of dark matter in the 1930s. Since then, obser- 
vational evidence has shown that dark matter behaves as a nonrelativistic fluid, \pdm\ *C e dm 
[3]. However, due to the lack of direct detection room is left for approaches where pd m ^ 
[3]. The dark matter equation of state can be inferred using gravitational lensing and kinematic 
techniques at astrophysical scales as proposed by Faber and Visser {5j . Relaxing the assumption 
Pdm = 0, dark matter in clusters of galaxies exhibits a negative equation of state parameter 
Wii m < [6j. For example, the Coma galaxy cluster, one of the four clusters analyzed in [6], 
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would be consistent with Wdm ~ —0.2. Similar constraints on the constant dark matter equation 
of state parameter can also be found at galactic [7j and cosmological scales [8]. 



If Wd m is negative enough to allow for the observed accelerated expansion of the Universe, 
one could abdicate dark energy. This is the idea behind unified dark datter (UDM) models. The 
UDM fluid guarantees a dark matter like behavior at high redshifts, where structure formation 
starts, and acts as a dark energy fluid at late times. The Chaplygin gas [9] and the viscous dark 
fluid (VDF) [2], with their "exotic" equation of state, realize this idea. Both models share similar 
background dynamics [10], whereas density perturbations of the Chaplygin gas are adiabatic (as 
CDM) and the VDF is intrinsically nonadiabatic [UJ. For this reason the Chaplygin gas mimics 
more properly cold dark matter, though it is phenomenologically motivated only. On the other 
hand, bulk viscosity is a well-known dissipative phenomenon in nature. Since real fluids are 
always subject to dissipation, it seems reasonable to assign this property to dark matter as well. 

The UDM models with bulk viscous pressure are competitive when confronted with super- 
novae la, baryon acoustic oscillations, age constraints and expansion rate measurements |12j and 
the matter power spectrum [TT]. However, the integrated Sachs- Wolfe (ISW) effect of such mod- 
els severely plagues models dominated by the VDF [131 EH E2] ■ The reason of such pathology 
is very simple: the bulk viscosity induces a large time variation of the gravitational potential at 
late times. 

In this work, we propose a viscous model that remains within the ACDM conception, where 
A drives the accelerated expansion of the Universe. We assume that CDM behaves as a real fluid 
equipped with bulk viscosity, viscous cold dark matter (vCDM). Within this approach we expect 
that the viscous effects do not cause a huge ISW effect, because (i) today's vCDM fractional 
energy density O v o ~ 0.25 (for the UDM model f2 v o ~ 1) and (ii) the equation of state of the 
viscous component assumes small negative values (in the UDM scenario we have w v o ~ — 1). 

In the next section we present the dynamics of the AvCDM model. Its background evolution 
is compared with current observational data in order to find the allowed viscosity of the vCDM 
component. The perturbative dynamics is also investigated. We calculate the ISW effect assum- 
ing the viscosity allowed by the background analysis. We also study the growth of subhorizon 
vCDM perturbations leading to the formation of dark matter halos during the matter dominated 
epoch. In some sense, we extend the findings concerning viscous UDM of our previous work 
[T5] . We discuss our results in the last section. 

2 Dynamics of the AvCDM model and observational constraints 
2.1 Background Dynamics 

The isotropic and homogeneous background dynamics of the AvCDM model resembles that 
of the flat ACDM model. The Hubble expansion rate H is given in terms of the fractional 
energy densities Oj, where the subscript i stands for baryons (b), radiaton (r), vCDM (v) and 
cosmological constant (A), 



In order to obtain the function O v (z) we have to specify the pressure of the vCDM and solve 
its continuity equation. 

Within relativistic thermodynamics the comoving frame is related to either the energy trans- 
port (Landau frame [16]) or to the particle number transport (Eckart frame [TT]). Here, we adopt 
the latter approach. For the homogeneous and isotropic background, dissipative effects manifest 
itself only through bulk viscosity. Shear viscosity and heat conduction could be relevant effects 
at the perturbative level, but we do not include them in this study. We assume that the vCDM 
fluid has a vanishing kinetic pressure and an intrinsic bulk viscous pressure 



h 2 = Hi [o b0 (i + zf + n r0 (i + zf + n v (z) + n A ] • 



(3) 



Pv = - 



(4) 
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The choice of £ seems to be the crucial aspect of any viscous model. In this work, we follow the 
recent results obtained in [131 02] an d stick to the same choice for the bulk viscous coefficient 



where £o an d v are constants and e v o is the density of the vCDM fluid today. This means that 
the current viscosity of such dark fluid is given by the parameter £q- 

For viscous UDM models, the ansatz ([5]) can lead to a large amplification of the ISW signal 
|13|. [T4l I15| . Fixing v = or v = —1/2, the ISW effect problem of these viscous cosmologies 
is less severe |15j . The quoted values for v have an explicit physical interpretation: the former 
means a constant bulk viscosity and the latter implies, for the one-fluid approximation (i.e. the 
bulk viscous fluid corresponds to the total density), the same background dynamics as seen in 
the ACDM model. We study the cases v = and v = — 1/2 and refer to them as models A and 
B, respectively. 

For model A, the continuity equation for the vCDM fluid reads 

(i + z)^^ - 3^v (z) + i [n l0 (i + zf + o b0 (i + zf + n v ( z ) + n A ] 1/2 = 0, (6) 



and for model B it becomes 

(i + z)^j^ - 3n v ( z ) + &li 2 n v (z)- 1/2 [o r0 (i + zf + n M (i + zf + n v (z) + n A ] 1/2 = o, 

(7) 

where the definition 

f-»gk ( 8) 
is valid for both models. The vCDM equation of state parameter w v = p v /ev today is 

Note that the CDM model is recovered if £ = 0. Fixing the values f2bo = 0.043 and f2 r o = 
8.32 x 10 -5 as given by the WMAP seven-year data pE], the remaing free parameters of our 
viscous models are £ and £l\ (with Q v q = 1 — fibo — ^ro — ^a)- 

2.2 Metric and density perturbations 

We assume a homogeneous, isotropic and spatially flat Universe, where the line element for 
scalar perturbations in the Newtonian gauge without anisotropic stress is 

ds 2 = a 2 (??) [- (1 + 2ip) dr] 2 + (1 - 2ij>) Sijcbtdx?] . (10) 

For convenience we introduce % = a' /a, where prime means derivative w.r.t. the conformal 
time rj. Let us first calculate the ISW effect [a net change in the energy of cosmic microwave 
background (CMB) photons produced by evolving potentials wells] of our viscous models. Given 
the gravitational potential ip it can be calculated by 

/AT\ P 

[=-) =2/ ifdq, (ii) 

V 1 J ISW Jr] r 

where the integration is performed along the photon trajectory from the conformal time at the 
recombination (r] r ) to the conformal time today (%)• 
The momentum constraint reads 

o tt2 2 

- k 2 i> - 3Hip' - 3H 2 ip = {n h A h + Q V A V } , (12) 
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where, for each component, we have defined the energy density contrast A = 8e/e. Since we are 
interested in the late time ISW effect, we have neglected the perturbations in the radiation fluid 
(Q r A r = 0). 

Additionally to (I12p we have the remaining perturbed Einstein equations 



k (V>' + Hip) 



¥*I«L {n h e h + (i + Wv )n v v } , 



(13) 



Hr)Q v 



in 
orL 



(14) 



where for each fluid we have defined the scalar velocity perturbation 9 by means of 8u,j = —kO ja 
with k being the wave number. 

Since there is no interaction between the different components of our model, each fluid obeys 
separately the perturbed energy-momentum conservation equations, 5T£,^ = (see [15] for a 
detailed derivation), which are used to close the set of equations. 

This set of equations must be solved numerically. Once we obtain the function ip, the ISW 
effect can be computed using ([11]) . 

Cosmological perturbation theory also allows us to investigate the formation of dark matter 
halos that attract baryons in order to form galaxies. The standard cosmological scenario gives 
rise to a hierarchical formation process where small structures form first. vCDM behaves as 
a pressureless fluid at the beginning of the matter dominated phase, i.e., w v (z ^> 0) ~ 0. 
However, the equation of state evolves in time departing from the pressureless behavior as the 
Universe expands. Additionally, there are also nonadiabatic contributions to the dynamics of 
density fluctuations. Assuming ([5]), the subhorizon vCDM perturbations obey the following 
Meszaros-like equation (see [15] for details): 



d 2 A, 
da 2 



+ 



a dH 
H~da 



+ 3 + A{a) + B{a)k 2 



dA x 
da 



+ 



+C(a) + D(a)k 2 



A v = P(a), (15) 



A(a) 



-6Wy + 



duu v 



2a dw v 



+ 



3u> v 



1 + w v da 1 + 2w v da 2(1 + w 



B(a) 



3a 2 H 2 (l + w v 



C(a) 



3Wy 



2(1 + w v ) 



3w v — 9w v 



3w 2 



1 + 



1 ^ a dH \ 

+ H~da) 



3a 



1 + 2w v \ dw 



1 + 



da 



+ 



6aw v dw v 
1 + 2w v da 



D(a) 



a 2 H 2 (l + w v 



. . dA v 

P{a) = —3vw Y a— h 3^w v A^ 

da 



1 9w v 
\ 

2 2 



+ 



-l-4w v + 2w 2 



dw v 



k 2 (l 



w v (l + w v )(l + 2w v ) da 3H 2 a 2 (l + w v ) 



where the scale factor was used as the dynamical parameter. The function P(a) contains the 
contributions from the perturbation of ^. Thus, for model A (v = 0) we have P(a) = 0. A 
remarkable point here is that the evolution of subhorizon vCDM perturbations is scale depen- 
dent. As shown in |15j . a viscosity large enough to accelerate the current Universe leads to a 
substantial suppression of growth at small scales. 
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Figure 1: Observational constraints on the parameter space (£, fl\) for model A (left) and 
model B (right). Solid lines are the contours of 2a confidence level. Long dashed lines are the 
age constraints (13 Gyr and 14 Gyr). Short dashed lines, from bottom to top, correspond to 
Q = and Q = 40% (see text for explanation). The horizontal dashed line sets the maximum 
allowed viscosity at 2a. 



2.3 Comparison with observations 

In Fig. 1 we show constraints on the free parameters £ and £1\ from different observational data 
sets. We constrain the background dynamics of our viscous models by means of the following 
observations: supernovae (SN) data (here we have used the Constitution sample [H]), the baryon 
acoustic oscillation (BAO) parameter A(z) |20j from the WiggleZ Dark Energy Survey [21] , and 
the position of the observed CMB peak l\ obtained by the WMAP project [22j that is related to 
the angular scale I a |23j . The solid lines in Fig. 1 are the 2a confidence level contours obtained 
from the likelihood function C oc exp(— \ 2 /2). The standard \ 2 statistics, 



N 



}th 7-)obs\2 



/ nth j~\c 

x 2 (l n A , Ho) = £ ( ' P ; , (is) 

^— < of 



measures the goodness of the fit. For each data set, with N data points, the theoretical value 
obtained within the AvCDM model D th is confronted with the observation D ohs . In order to 
obtain the bidimensional likelihood contours shown in Fig. 1, we marginalize Hq with a flat prior 
[0 < i^o(Km _1 sMpc) < 100]. Long dashed lines are age constraints corresponding to 13 Gyr 
and 14 Gyr. The gray filled areas in both panels of Fig. 1 correspond to the "concordance" 
parameter values allowed at 2a confidence level. The best fit occurs at £ = 0, but it is possible 
to establish an upper bound (at 2a) to the viscosity parameter £. For the model A (B) this 
value is £ < 0.24 (0.31) as seen in the horizontal dashed line in the left (right) panel in Fig. 1. 

With the equations for the potential, derived in Sec. 2, we compute the ISW effect of both 
viscous models and compare it with the prediction of a flat ACDM model (the fiducial cosmology 
adopted here has parameters Hq = 72 km/s/Mpc and f2 m o = 0.266, as suggested by the WMAP 
seven-year analysis |18]). The short dashed lines in Fig. 1 correspond to relative amplifications 
(Q) of the ISW effect calculated as 



/ AT \ AvCDM 

[~T~J 

/ AT xACD\l 



O = v LMSW 1 (17 \ 

^~ / AT xACDM ' 
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If Q > (< 0) the AvCDM model produces more (less) temperature variation to the CMB 
photons via the ISW effect than the fiducial ACDM model. Such an analysis of the ISW effect 
has been proposed in [23] and used in [15]. As seen in Fig. 1, the parameters allowed by the 
background correspond to values between Q = 0% and Q = 40%. This is a reduced ISW effect, 
when compared with the viscous UDM models where the amplification could reach Q ~ 120% 

m- 

A successful structure formation process (within the hierarchical scenario) is achieved when 
small structures (formed from the smallest halos) merge into large matter agglomerations. At 
the moment of matter-radiation equality (z eq ~ 3000) typical CDM subhorizon perturbations 
grow oc a. Equation (|15p also has this solution in the limit w v = 0. However, the vCDM 
(w v / 0) perturbation evolution is scale dependent. 

In Fig. 2, we plot the evolution of the density contrast A for both viscous models for scales 
in the range k = lOOO/i Mpc -1 (dwarf galaxies) to k = 0.2h Mpc -1 (galaxy clusters). In 
each panel the solid line is the standard CDM growth oc a. The dashed lines correspond to 
the viscous models for different values of the viscosity coefficient. The initial conditions, i.e. 
the power spectrum at the matter-radiation equality, are set using the CAMB code [25] . This 
provides the correct amplitude for each k mode at z eq and helps us to identify the onset of the 
nonlinear structure formation A(z n i) = 1. Considering £ ~ 0.2, as obtained above, the viscous 
effects would suppress the halo growth well before z e q. Viscous dark halos at cluster scales 
(k = 0.2/iMpc -1 ) are able to follow the typical CDM growth only if £ < 10~ 6 . Smaller scales 
place stronger constraints on the bulk viscosity, for dwarf galaxies £ < 10~ n in model A and 
| < 10~ 10 in model B. 

In Fig. 3 we implement a more conservative analysis. It shows the maximum viscosity 
allowed in order to reproduce the standard CDM perturbation growth until A = 1. The left 
panel displays the constraints on the dimensionless quantity £. The right panel shows the same 
results in SI units. The interpretation of these plots is as follows: For the wave numbers shown 
on the horizontal axis, linear perturbations of the vCDM fluid are identical to standard CDM 
if the viscosity is below the indicated critical values. For larger values, structure formation is 
affected. 

3 Conclusions 

In order to investigate the dissipative properties of dark matter, we have proposed the AvCDM 
model, where cold dark matter is modeled as a dissipative fluid equipped with bulk viscous 
(negative) pressure (the vCDM fluid). Differently from UDM viscous models, where the viscous 
dark fluid is responsible for the accelerated expansion of the Universe, here we allow the existence 
of a nonvanishing cosmological constant. This means that the vCDM equation of state can be 
only slightly negative. This possibility is motivated by recent studies of galaxy clusters [6]. 

The free parameters of our model are f^A and the dimensionless viscosity parameter £;. The 
limiting cases f^A = and £ = correspond to the viscous UDM model (the VDF model) and 
the ACDM model, respectively. The comparison with observations (see Fig. 1) has revealed 
that the unified model is strongly disfavored. The contour (at 2a) for the WiggleZ BAO data 
together with the position of the CMB first peak is decisive to put an upper limit on the dark 
matter viscosity £. For both vCDM dynamics studied here, namely models A {y = 0) and B 
(y = —1/2), we obtain £ < 0.2 (at 2a) corresponding to £o < 10 Pa-s in SI units. 

For the sake of comparison, the bulk viscosity coefficient of water at atmospheric pressure 
and at room temperature at 25 °C is 2.5 x 10~ 3 Pa-s [27] : thus we see that constraints from the 
dynamics of the homogeneous and isotropic background are rather poor. 

For the allowed (at 2a) model parameters, the AvCDM does not produce a large amplification 
of the integrated Sachs- Wolfe signal as in the viscous UDM model case. In the latter approach 
the background preferred data produces Q > 120% [15]. The current estimations of the error 
bars in the cross correlation of galaxy density and CMB temperature (used to measure the ISW 
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Figure 2: Growth of subhorizon density perturbations for different scales in the range k = lOOO/i 
Mpc -1 (dwarf galaxy), in the bottom panels, to k = 0.2h Mpc -1 (galaxy clusters) in the top 
panels. Left (right) panels correspond to model A (B). The solid line corresponds to the standard 
CDM A oc a result. Dashed lines correspond to the viscous CDM growth for various values of 
the parameters £ as indicated in the panels. 7 
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Figure 3: Maximum viscosity allowed following the requirement A v cdm(^) = Acdm(-z) for 
z > Zni- For £ > £, m ax and £o > (omox! structure formation is affected by viscous phenomena. 



effect) are still compatible with models where Q ~ 100% [26J. This means that amplifications 
of the order Q < 40%, as found here for the AvCDM model, are compatible with current 
observations. Radio surveys in the near future will reduce the error by a factor of 5 and thus 
should be able to improve the limit to the Q ~ 20% level (see figure 9 of |28j). 

The AvCDM is much more tightly constrained by the analysis of the growth of subhorizon 
density perturbations. The production of small halos is a fundamental aspect of the hierarchical 
structure formation process. When even a small bulk viscous (negative) pressure is present, the 
growth of the smallest dark matter halos is suppressed as shown in Fig. 2. The suppression 
in model A (a constant bulk viscosity) is less severe than model B (which during the matter 
dominated phase is equivalent to a constant negative pressure). As seen in Fig. 2, vCDM is 
able to form galactic dark halos, within the hierarchical scenario, only if £ <C 0.2. For dwarf 
galaxy scales (~ 1 kpc) the growth of the density contrast is similar to the standard CDM if the 
viscosity values are reduced to £ < 10 -11 . However, note that avoiding the formation of dwarf 
galaxies could provide a solution to the problem of the missing satellites in the standard CDM 
scenario [29J. Assuming that we have to guarantee, at least, the existence of 10 kpc galaxies, 
the viscosity parameter is set to £ < 10 -9 or £o < 10 _1 Pa-s in SI units. This result is even 
smaller than the limit £ < 7.38 x 10 3 Pa-s, obtained in [30], where the authors study the Jeans 
mechanism for the bulk viscous fluid at the recombination era. 

Since the parameter £ is of the same order of magnitude as the equation of state parameter 
today w v o [see ©], the results mentioned in the introduction, i.e w^m ~ —0.2 for the Coma 
cluster, are challenged by our analysis. 

Our conclusions are limited by the linear analysis here performed. This formalism breaks 
down when A ~ 1 and numerical simulations would be required to predict the final clustering 
patterns. The study of the nonlinear collapse of a vCDM component can be decisive to clarify 
whether virialized viscous structures survive at the end of the matter dominated epoch. A 
further limitation of our study relies on the assumption that the coefficient £ obeys the ansatz 
([5]). A microscopic model for the bulk viscosity, as one finds in the "dark goo" model [31], would 
yield to a more appropriate coefficient £. 

To summarize, we have shown that the strongest constraints on the dissipation of dark matter 
come from the study of hierarchical structure formation at the smallest observable scales. As 
dwarf galaxies are observed, their existence implies that £o < 10 -3 Pa-s, a value as large as the 
bulk viscosity of water at normal conditions on Earth. 
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